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CALCULATION  01?  THE  ATTENUATION  COEFFICIENTS  OF  0.6-14  pm  WAVES 

PASSING  THROUGH  FOG 


I.  L.  Zel'manovioh,  L.  M.  Lobkova,  and 
Ye.  R.  Milyutin 

Statement  of  the  Prob lem 


The  aerosol  (drops  or  bits  of  ice  in  precipitation  and  fog,  smoke, 
and  dust  particles)  is  the  basic  variable  part  of  the  troposphere 
determining  its  optical  properties  in  the  visible  and  infrared  areas  of 
the  spectrum.  The  composition,  shape,  and  structure  of  the  aerosol  are 
the  factors  that  determine  the  transparency  of  the  atmosphere. 

On  the  basis  of  physical  presentations,  the  attenuation  coefficient 
of  electromagnetic-  waves ,  a  1/unit  of  length,  can  be  presented  as 

«  •  «.  +  °a  +  “p  *  “aer-  (1) 

Here  o,  is  the  coefficient  of  molecular  scattering;  a  is  the 

5  a 

coefficient  of  molecular  absorption;  ap  is  the  attenuation  coefficient 
caused  by  the  radiation  pressure  (due  to  motion  of  drops  or  particles); 
aaer  is  the  coefficient  aerosol  attenuation. 

The  coefficient  of  the  molecular  scattering  for  the  visible  portion 
of  the  spectrum  is  negligibly  small.  For  example,  for  the  wavelength 
A  *  5500  A  a  s  0.012  1/km  [1]. 


FTD-HT-2 3-59 5-70 


1 


Molecular  absorption  in  the  windows  of  atmospheric  transparency 
can  also  be  disregarded  fl]. 

The  absorption  and  dispersion  of  waves  of  the  optical  hand 
connected  with  the  radiation  pressure  must  be  considered  when  the  force 
of  the  pressure  is  sufficient  to  bring  the  water  drops  or  particles  in 
the  atmosphere  into  motion.  In  their  movement  through  air  the  drops 
(particles)  generate  friction,  which  is  another  cause  for  the  attenuation 
of  the  field  and  the  oscillation  of  the  dipole  moment  in  a  unit  volume 
of  air  caused  by  the  motion  of  the  drops  will  be  another  cause  of 
scattering.  According  to  work  [2],  the  radiation  pressure  effect  occurs 
only  with  very  high  radiating  energies  and  small  drops  (r  >  1  pm). 

Thus,  the  effect  of  radiation  pressure  on  the  attenuation  of  the 
wave  can  be  disregarded. 

Consequently,  the  attenuation  of  waves  of  the  optical  range  will 

be  determined  basically  by  the  presence  ef  aerosol  in  the  atmosphere. 

Prom  now  on  we  will  assume  that  a  *  <x  . 

&er 

This  work  is  concerned  with  the  calculation  of  the  attenuation 
coefficient  In  mist  and  fog,  using  exact  and  approximate  formulas  of 
the  scattering  theory,  and  also  with  the  comparison  of  the  obtained 
results. 

Determination  of  the  Attenuation  Coefficient  in  Mist  and  Fog 
Using  Exact  and  Approximate  Formulas 

The  attenuation  coefficient  of  the  electromagnetic  oscillations, 
a  dB/km,  scattered  in  mist  and  fog  can  be  calculated  using  .the  formula 
[31 

a  —  4,34  •  101  i'c(/u,  -,)N{r)dr.  ( 2 ) 

0 

Here  N(r)  is  the  function  cf  distribution  of  the  drops  with  respect  to 
size;  N(r)dr  is  the  relative  number  of  particles  whose  size  is  within 
the  limits  r,  +  dr;  a(m  ar.d  p)  is  the  coefficient  of  attenuation  on 
a  spherical  drop  (unit  of  area);  m  =  n-.i<  is  the  relative  complex  . 
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index  of  refraction  on  the  water-air  boundary;  n  is  the  index  of 
refraction;  k  is  the  absorption  coefficient;  p  =  is  the  wave 
parameter;  r  is  the  radius  of  the  drop;  A  is  the  wavelength  of  the 
Incident  radiation. 

Let  us  examine  the  values  a(m,  p)  and  N(r)  which  enter  into 
expression  (2).  Experimentally  it  was  established  that  the  distribution 
curve  N(r)  of  drops  with  respect  to  size  in  mist  and  fog  is  described 
well  by  the  four-parameter  family  of  functions  [4  and  5] 

N(r)  =  Ar+e->'\  (3) 

where  A  is  the  normalization  parameter;  u,  $,  and  y  are  curve  parameters. 

The  four-parameter  family  of  curves  (3)  is  very  general  and  permits 
us  to  describe  with  the  necessary  accuracy  practically  any  single-apex 
distribution  curve  and,  with  summation  of  the  expressions  of  type  (3), 
also  any  multi-apex  curve  of  distribution.  Giving  determined  values  to 
quantities  y  and  we  obtain  the  known  distributions  of  the  aerosols. 
For  example,  in  1930  I.  Rokar  proposed  for  the  aerosol  the  distribution 
(3)  with  y  *  2  and  y  *  1: 

jV{r)*=Ar-e-t'.  (*0 

In  1952  H.  Young  found  that  in  the  surface  layer  £4] 

N(r)  =  (5) 

r 

in  this  case  we  have  8  a  0  and  y  **  -4. 

If  we  standarize  distribution  (4)  to  N  —  the  number  of  particles 
in  1  cm^  -  we  will  obtain 


A 


(6) 
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Q.c;o^.7n 

-  »•*  *- -/  ✓✓  ✓  i  - 


Then 


:a-i 


(7) 


u  +  l 


Here  is  the  gamma-function  from  the  argument  — - —  and  parameter 

8  is  determined  by  the  expression 


* 

where  r  is  the  average  radius  of  the  drop. 


(8) 


The  majority  of  the  experiments  indicate  that  we  can  assume  y  *  1. 
Parameters  y  and  r  are  given  to  describe  the  .microstructures  of  the 
mist  and  fog. 


The  calculation  of  the  coefficients  of  attenuation  on  a  spherical 
drop  a(m,  p)  can  be  accomplished  using  the  exact  and  approximate 
formulas  for  scattering  of  light  [3  and  6].  According  to  Mie*s  theory, 
the  coefficient  of  attenuation  on  a  spherical  drop  in  air  (n  *  1  and 
k  *  0)  is  determined  by  expression  [6] 

*  (».  fO  -  £  Im  2J  /  (/+1K-1)1  (c,-bt).  (9 ) 

1 

Prom  (9)  it  is  evident  that  the  diffracted  field  is  presented  in 
the  form  of  the  sum  of  the  individual  partial  waves  (cz  —  are  electrical 
waves  and  b^  are  magnetic  waves).  The  intensity  of  excitation  of  the 
Ith  partial  wave  is  determined . by  values  c^  and  b^  which  depend 
essentially  on  p  and  on  m. 


In  turn,  the  partial  amplitudes  c^  and  b^  can  be  determined  by  the 
relationships 


ct-i  (-D' 


it  -  i 

m-r  1) 


?)  -  « >!'/(«  ?V 


*|««/(- 1)‘J‘ 


2  /  -l  i 

7<v~r 


?)-w  ^(?)a-t\'i  p) 
5;'(p)'r;  («  ?)— «  «/  (f)yl(m  p) 


(10) 
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f 

s* 

4 

f 


/ 


Coefficients  and  bj  are  determined  by  the  cylindrical  Bessel 
and.  Hankel  functions  with  half-integer  indexes.  These  functions  are 
determined  by  the  equations 


=  |/"  y-  i  (z)  —  ^ y  • 


(11) 


To  determine  the  function  ^(z),  one  of  the  relationships  between 
the  solutions  of  the  Bessel  equation  is  used: 


(z)“<Mz)+*Xj(2), 


(12) 


where 


7.i  (z)  “*  (—  *=  -*]/*• y  Y(t+'i,)(z)‘ 


(13) 


To  calculate  functions  ^j(z)  and  x^(z)  30(1  their  derivatives,  we 
use  the  recurrence  relationships  of  the  type 


»*»(*) 


(14) 

<(*).- W«) — f^z).  (15)^ 

s' 

When  relationships  (14)  and  (lb)  are  used  the  initial  functions  are 


♦#{z)«sinar;  •{»,(*)  - 


kin  t 


-cos  z; 


Xo(z)-cosz;  &(*)<■«•  sin 


(16) 


However,  the  calculation  of  functions  ♦j(p),  ^(pm),  and  £j(p)  and 
their  derivatives  directly  by  recurrence  relationships  (14)  and  (15) 
leads  to  loss  of  accuracy  for  the  large  values  of  l.  For  example,  if 
the  values  of  functions  $q(z),  ^(z) ,  Xq(z)»  and  X^(z)  are  known  with 
an  accuracy  to  10”®,  then  ^7(z)  and  57(z;  with  k  >0.1  and  l  *  80-100 

^  -1  o 

will  have  an  accuracy  up  to  10  -10  .  The  given  problem  was  subjected 
to  a  detailed  3tudy.  It  turned  out  that  for  the  direct  calculation  of 
function  tj^(rap)  (p  >  50 j  n  <  1.16;  <  >  1.16)  with  the  aid  of  summing 
of  the  series  [7] 
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(17) 


.  ,  .  »/sW?  .  /  \  ‘  fr'  ( '  m  V 

?;(«;>)=  /  •/.(«»*)  —  /  — r-  _  JT 

f  '  -  .S  0  °  ’ 


(-1)V  /  MlV 

6'! i‘(/  —  A'  ••  3/.)\  2  ' 


a  great  deal  of  machine  time  would  be  required. 


Consequently,  the  method  of  direct  calculation  of  function  ^..(mp) 

fs 

using  expression  (17)  is  inappropriate. 


For  calculation  of  the  function  ^(mp),  the  so-called  method  of 
dispersion  [8  and  9]  was  used;  this  ensures  the  calculation  of  the 
attenuation  coefficient  at  large  p  and  k  >  0.1  to  four  places.  Such 
accuracy  fully  satisfies  all  practical  questions. 


In  summing  expression  (9),  the  value  of  S  was  chosen  thus:  for 
p  =■  0.5-1. 5,  S  -  3;  For  p  =  2-4,  S  =  5-7;  for  p  =  8-12,  S  *  14-18;  and 
for  p  >  13,  S  =  (l.5-1.2)p. 


Calculation  of  the  attenuation  coefficient  a  is  connected,  as  is 
indicated  above,  with  great  difficulties;  moreover,  calculations  using 
formula  (2)  can  only  be  accomplished  by  the  numerical  integration 
method. 


The  calculation  of  attenuation  coefficients  can  be  simplified 
considerably  by  using  tables  of  light  scattei’ing  [10].  However,  the  use 
of  these  tables  does  not  eliminate  the  problem  of  using  approximation 
formulas  to  calculate  attenuation  coefficients  a. 

The  coefficient  of  attenuation  on  a  drop  of  water,  o(n,  p),  in  the 
absence  of  absorption  (x  *=  0)  is  determined  by  [3] 

o(/f,  p)«srJ  |2~  -cos*)].  (18) 

Taking  into  account  the  absorption  (x  ^  0), 

a(m,  p)  =  -:r3|2  —  4e-*w,-~j~sin(.x  —  v)  — 

— 4e— rtc  [-"'-J’cos  {x — 2^)  4-  4|  —-y- )  cos2*j.  ( 19  ) 

x  —  2,o  {«—!),  tg  v  =  ~-j ■{ - . 
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Formulas  (18)  and  (19)  are  valid  with  the  following  values: 
n  =  1.05-2  and  (m-1)  <  1. 

Substituting  in  expression  (2)  functions  (7)  and  (18)  or  (19),  we 
obtain  the  formula  for  the  determination  of  attenuation  coefficient 
a  1/km  both  in  a  nonabsorbing  poly dispersed  sol 


*  PN  •  10»  4,  / 

* — Rim rg/1' 


2r(i*  4-  3) 


/, «  -2. . sin[(i‘  +  2)  arc  tg  f] , 

IP  +  (2c>-i  2 

i  roi  +  n' 


i 


/»* 

r(|i. 


y»+l 


Ci  £+‘ 

IP  +  (2e):)  2 

c  —  1), 


“r 'cos  [(n  + 1)  arc  tg  y]. 


(20) 


and  in  an  absorbing  polydispersed  sol 


/, 


/<■=  - 
/»- 


a 


- 10* 

~W+rr 


l  '  ■ 

' 

,  _  Kh-tgg.  3in[(i*+  2)  arc  tg-^-j, 

(*»  +  9s)  * 

47»;»ln’  B£±|L.  cos  [fa + 2)  arc  tg  , 

fa’  +  ff5" 

Aco^Fcoejv.  cos  [((a  4*  Darctgi-). 

(*  r  'T')1*’ 

ifa  + +L.Sin[fa+  Dsrctgy], 

-j-  95) 

,  «cos  >)=  co*  2»  Ifa-t- 1)_ 

<o  • 

6  «=  -  11  /**  tg v,  «r /  +  8. 
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Comparison  of  the  Attenuation  Coefficients  in  Mist  and  Fog 
Calculated  by  Exact  and  Approximate  formulas 
of  Light  Diffusion 

Comparative  calculations  (see  Table  2)  were  made  for  several 
operating  waves  of  a  laser  which  lie  in  the  spectral  windows  of 
atmospheric  transparency  [11]. 

The  optical  properties  of  liquid  water  for  the  selected  laser 
emission  waves  were  obtained  by  the  graphic  interpolation  of  data 
presented  in  works  [6,  12].  These  works  give  values  of  n  and  k  up  to 
X  <  9  pm  at  t  =  20°C  and  for  X  >  9  pm  at  t  =  l8°C. 

On  the  basis  of  the  experimental  data  on  the  microstructure  of  the 
mist  and  fog  the  following  values  for  the  radii  of  drops  were  selected: 
0.2  (0.2)  1  (0.4)  5  (0.5)  14  (2)  60  (number  in  parentheses  —  step  of 
calculation),  and  the  mean  drop  radii  were:  1  (mist)  2,  3t  5,  7,  10, 
and  15  pm.  The  given  drop  sizes  basically  characterise  the 
microstructure  of  dense  mist  and  fog. 

The  attenuation  coefficients  for  the  mist  and  fog  per  1/km  were 
determined  with  the  condition  that  the  concentration  of  drops  N  »  1  in 

3 

1  cm  of  air. 

To  determine  the  effect  of  the  parameter  p  on  the  attenuation 
coefficient,  a  numerical  experiment  was  carried  out  on  the  wavelength 
X  «  O.63  pm  of  emission  of  a  He-Ne  laser;  the  optical  constants  of  the 
water  —  air  medium  were  n  »  1.3318;  <  ■  0. 

The  calculations  were  done  by  the  formula 


« =>=  101  V  c(/«,  5 >}N(r)  A  r,  ( 22 ) 

taking  into  account  expression  (9)  and  distributions  (7)  for  p  equaling 
0,  2,  4,  6,  8,  10,  and  20. 

The  tabulated  values  of  a  for  1  km  are  presented  in  Table  1.. 
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Table  1 


X--6328  A;  /i- 1,3318; 

? 

ym 

i 

JAltlO  j 

*=2 

;a  =  4  j 

i 

i1— G 

1 

I 

,=8  j 
1  *  I 

1 

j  ;ji~10  ' 

;i=-20 

1 

1,405  -2  , 

9,700  -3 

i 

8,744  -3  | 

8,257  -3- 

7,923  -3 

|  7,661  -3 

6,811  -5 

2 

5,359  -2 

3,679  —2 

3,335  -2 

3,188  -2 

3,108  -2 

3,056  -2 

2,936  —2 

a 

1,168  -! 

8,101  -2 

7,341  -3 

7,018  -2 

6.842  -2 

6,732  -2 

6,508  -2 

5 

3,164  ~1 

2,168  -1 

1,985  —1 

1,901  -1 

1,851  -1 

1,819  -1 

1,746  -1 

7 

6,166  -1 

4,135  -1 

3,754  —1 

3,620  —1 

3,551  -1 

3,306  — 1 

■  3,397  — : 

10 

1,204 

8,324  -1 

7,413  —1 

7,033  -1 

6,837  -1 

6.726  -1 

I  6.538  --1 

15 

2,235 

1.893 

1,696 

1,601 

1.544 

1,505 

1,406 

Note.  Here  and  in  Table  2,  the  fifth  number, 
separated  from  the  preceding  oy  spacing,  designates 
the  coefficient  on  base  ten  by  which  the  four 
digital  number  preceding  it  will  be  multiplied. 

•Hie  variation  of  the  attenuation  coefficient  in  fogs  as  a  function 
of  parameter  y  *  2-20  does  not  exceed  20%,  which  indicates  that  the 
effect  of  the  parameter  y  on  the  quantity  a  is  small.  If  on  the  basis 
of  [33  we  assume  that  y  •  6-8,  then  in  this  case  with  variation  of  y 
from  2  to  8,  a  will  not  vary  more  than  10-15%. 

In  calculating  attenuation  coefficients  in  mist  and  fog  it  was 
assumed  that  y  2  and  r  equaled  ym  1,  2,  3,  5,  7,  10,  and  15. 

The  use  of  distribution  with  y  »  2  is  Justified  by  our  effort  to 
obtain  a  certain  reserve  coefficient  during  the  determination  <*£  the 
range  of  action  of-  optical  means  of  communication  and  location. 

Table  2  shows  the  values  of  the  poly dispersing  attenuation 
coefficients  calculated  by  the  approximate  formulas  (20)  and  (21) 
with  N  «  1  per  1  cm^  of  air. 

The  data  of  Table  2  shows  that  the  error  in  the  determination  of 
the  attenuation  coefficients  using  approximate  formulas  does  not  exceed 
~l8%  in  comparison  with  the  calculation  results  obtained  by  the  exact 
formulas i  this  is  entirely  acceptable  for  carrying  out  engineering 
calculations . 
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Table  2 


1 

X  inn  j 

■  1 

?  MKM  ! 

1 

■l  if  KM 

apr* 

*  i 

Relative 

error 

X 

0.6328  ! 

1,3318  j 

> 

0 

1  : 

9.700  —3  ! 

9.203  —3  i 

5.1 

! 

_  f 

3,679  -2  j 

3.590  -2  ! 

2,7 

i 

! 

i 

3,101  -2  j 

7,820  —2  ' 

3.4 

1 

t 

„  1 

J  i 

2,162  -1  j 

2,199  —2  ! 

5.0 

I 

j  * 

7  ' 

4.135  ~1  ; 

4.341  -1 

6.4 

i 

!0  ; 

'  8,324  -1  j 

8,855  -1  j 

5.5 

15  i 

1,993  j 

1,997 

1.4 

1.0621 

1,3247  | 

2.8- 10"*  i 

‘  | 

1.063  —2 

1,024  -2 

3.7 

i 

i  ! 

2  , 

3.853  -2  1 

3,430  -2 

11.0 

* 

»  ! 

t 

3 

8560  —2  | 

7520  -2 

10.0 

5*  ; 

2,224  ■*!  j 

2,025  -1 

9,0 

7  j 

4,247  -1  j 

3,994  —1 

5.4 

1 

1 

i 

10  ! 

8,693  -1  ! 

8,218  -1 

5,5 

! 

1 

i  is  | 

1,999  j 

1,852  1 

7.4 

2,095 U  i 

1,302 

|  155- 10"3  | 
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14.3 
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3 
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1 

• 

5 
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63 

1 

7 
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l 

\ 
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13 

i 

15 
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14 
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1 
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1 

; 

2 
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3395  -2 

183  . 

3 

1,060-1 

9.417  —2 

U3 

! 

s 

254$  -1 

2315  -1 

93 

' 

•f 

•  * 

4423  -J. 

4.450-1 

33 

*  ! 
! 

. 

10 
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13 

» 
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!  14 
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853 -.10-1 
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.  1 

! .  ] 
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*> 
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10 
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•  i 

! 
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1  2.034 

!  1,997 

i  .  12 

Note.  The  asterisk  denotes  the  r  values  after 
which  the  calculations  of  o  were  carried  cut  with 
less  accuracy. 
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Figure  1  shows  the  ratio  of  the  attenuation  coefficients  of  various 
wavelengths  in  mist  and  fog  to  the  attenuation  coefficient  for  the 
wavelength  X  *  0.6328  yni. 


et 


Fig.  1 


The  direct  measurement  of  the  quantity  of  drops  in  mist  and  fog  is 
a  tedious  operation  and  does  not  permit  fast  calculation  of  the 
attenuation  coefficient.  At  the  same  time,  with  the  aid  of  radars  the 
water  content  of  fog  located  any  place  along  the  route  can  be  determined 
from  a  distance.  For  this  reason  in  work  [13]  the  relationships 
connecting  the  attenuation  coefficient  with  the  water  content  of  the 
fog  were  obtained. 


Conclusions 


1.  The  results  of  calculations  using  the  approximate  and  exact 
formulas  coincide  with  an  accuracy  of  up  to  20$. 

2.  In  the  range  0.6-14  yra,  electromagnetic  waves  in  the  10-12  ym 
region  pass  through  fog  with  the  least  attenuation. 

3.  In  heavy,  optically  dense  fogs  (for  short  X)  the  coefficient 
of  attenuation  is  not  selective. 
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